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ABSTRACT 

In X-ray binaries, rapid variability in X-ray flux of greater than an order of mag- 
nitude on time-scales of a day or less appears to be a signature of wind accretion from 
a supergiant companion. When the variability takes the form of rare, brief, bright out- 
bursts with only faint emission between them, the systems are called Supergiant Fast 
X-ray Transients (SFXTs). We present data from twice- weekly scans of the Galactic 
bulge by the Rossi X-ray Timing Explorer (RXTE) that allow us to compare the be- 
haviour of known SFXTs and possible SFXT candidates with the persistently bright 
supergiant X-ray binary 4U 1700-377. We independently confirm the orbital periods 
reported by other groups for SFXTs SAX J1818.6-1703 and IGR J17544-2619. The 
new data do not independently reproduce the orbital period reported for XTE J1739- 
302, but slightly improve the significance of the original result when the data are 
combined. The bulge source XTE J1743-363 shows a combination of fast variability 
and a long-term decline in activity, the latter behaviour not being characteristic of 
supergiant X-ray binaries. A far-red spectrum of the companion suggests that it is a 
symbiotic neutron star binary rather than a high-mass binary, and the reddest known 
of this class: the spectral type is approximately M8 III. 

Key words: X-rays: binaries - stars: supergiants - accretion, accretion discs - stars: 
neutron - binaries: symbiotic 



1 INTRODUCTION 

Blue supergiant stars were first identified with fast X- 
ray transients d iscove r ed by AS C A (A X 1845.0-0433, 



lYamauchi et al.l (|l995l ): [Coe et al.l (Il996l) ) and later b y 
RXTE (XTE J1739-302. ISmith et all (|l998l ): ISmithl (|2003h ). 
As more began to be disco vered with INTEGRAL , they 
were recognised as a class (ISmithI |2004 lin't Zandl l2005l : 
INegueruela et al.l l2006l : ISguera et al.l (2005) and named Su- 
pergiant Fast X-ray Transients (SFXTs) (Nc gueruela et ah! 
l2006l ). 

SFXTs are thought to be a subclass of X-ray binary 
in which a neutron star accretes the wind of a supergiant 
companion. Some binaries with these two components have 
been known since the 1970s as very bright and variable X-ray 
sources (e.g.. Vela X-1, 4U 1700-377), but the SFXTs spend 
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most of their time at low luminosities of 10"^^ to 10^^ erg s~^, 
with only very brief excursions to outburst luminosities of 
up to a few times 10'^'' erg s^^. These outbursts generally 
last for seve r al hou r s, with subpeaks last ing for minutes (e.g. 
ISidoh et all (|2009l ): iRampv et all (|2009l )). 

The cause of this behaviour is not known. Clumps 
in the wind of the companion can explain rapid variabil- 
ity w i thin outbursts llin'tZandll2005l : I Walter fc Zurita HerasI 



I2OO7I : INegueruela et al.ll2008l ) but are less successful at ex 
plaining long periods of quiescence. In some ca ses a highly 
eccentric or wide orbit may be responsible jSidoli et al.l 
|2003), but some SFXTs have been shown to have orbits 
as short as the persistently bright sources (see Table 1). 
Rej ect ion of the wind material by the magnetic propeller 
effect (iGrebenev fc Sunvae ^ |2007I :[b ozzo. Falanga fc Stellal 
I2OO8I: iLi fc Zhand |201ll) has been suggested. Depending 
on the particular systems studied and the assumptions 
made about the accretion physics, consideration of the pro- 
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J ieller effect has led to both h igh lower limits (> fO^'^ G) 
Bozzo. Falanga fc Stella 20081 ) and low upper limits (> 
10" G) (|Li fc Zhanel \20n\ ) on the magnetic field. Unex- 
pected structure in the companion wind itself, e.g. very rare 
clumps with a very high density contrast, has not yet been 
ruled out. 

Some sources have been referred to as intermedi- 
ate between SFXTs and the familiar persistent super- 
giant/neutron star binaries, but this category has never been 
universally defined. The most common criterion for "inter- 
mediateness" mentioned in the literature is the dynamic 
range between the highest and lowest luminosities observed. 
"Intermediate" sources by this criterion have been defined 
as either those with a dynamic range of <100 but >20 or 
so (jWalter fc Zurita He ras 2007; Clark ct a l. 2010) or those 
with a range of 100-1000 (|Sguera et al.ll201ll : iRomano et al.1 
|2012| ). The maximum dynamic range of any source - both 
the SFXTs and the classical wind-accreting HMXBs such 
as Vel X-1 - can only increase with the number of observa- 
tions made. Thus it is not surprising that the most recent 
of these references give the higher range for intermediate 
systems, since the known dyn amic ranges o f all the objects 
have increased in recent years. IChatvl (|201G| ) adds two addi- 
tional criteria to the definition of intermediateness: a higher 
average luminosity than the canonical SFXTs, and a longer 
average duration for individual flares. [Romano et al.l (|2012l ) 
add the criterion that the lowest luminosity frequently ob- 
served should be close to 10^^ erg s~^ for intermediate sys- 
tems, and 1 to 2 orders of magnitude lower for true SFXTs. 
iRomano et all (l201ll ) have begun to make the SFXT vs. 
intermediate classification systematic for the first time by 
interpreting histograms of the luminosity of three sources, 
XTE J1739-302, IGR J17544-2619, and IGR J16479-4514 
made from monitoring data from the sensitive Swift X-Ray 
Telescope (XRT). The difference between true SFXT and 
intermediate systems in that work is based on the fraction 
of time spent at the lowest luminosities. 



2 DATA AND ANALYSIS 

We present data fro m the Galactic bulge scan programme 
l| Swank et al.l 1200 ll ) of the Proportional Counter Array 
(PCA) in strument on the R ossi X-ray Timing Explorer 
(RXTE) jjahoda et al.1 120061 '). PCA is a large-area X-ray 
detector with a collimated but non-imaging field of view of 
1° full width at half maximum. Twice a week, PCA scans 
the inner Galaxy in a raster pattern, recording the instan- 
taneous 2-10 keV count rate of a large catalogue of sources. 
The bulge scan data are publicly available. 

In this paper we show data from the two 'proto- 
type' SFXTs XTE J1739-302 and IGR J17544-2619 (Fig. 1 
& 2), the persistently bright wind accretor 4U 1700-377 
(Fig. 3), three supergiant binaries that are not as well stud- 
ied and have been considered either SFXTs or 'interme- 
diate' systems (SAX J1818.6-1703, AX J1845. 0-0433, and 
IGR J16479-4514, Fig. 4, 5, and 6, respectively), and finally 
two fast transients whose companion type and nature are not 
known: XTE J1743-363 and XTE J1901-F014 (Fig. 7 and 8). 



^ http:/ /Iheawww. gsfc.nasa.gov/users/craigm/galscan/main. html 



For each source we present four complementary views of the 
data: the 2-10 keV light curve from up to 12 years of bulge 
scans, a Lomb-normalised periodogram used to search for or- 
bital periodicity, an integral histogram of the luminosities, 
and the autocorrelation function of the light curve. 

The light curves can be considered reliable at relatively 
high count rates, but near zero counts the results are com- 
plicated by the systematic limitations of the technique of 
fitting the count-rate time profiles of the raster scans. For 
example, systematic uncertainties in the profile of the Galac- 
tic 'diffuse' emission (much of which may be an ensemble 
of many faint ca t aclysm ic variables and other stars - see 
iRevnivtsev et al.l Hooi)) translate to uncertainties in the 
lowest source count rate. This problem is particularly severe 
for AX J1845. 0-0433 (Fig. 5a), where the bifurcated ap- 
pearance of the baseline values is caused by the alternation 
of scans parallel and perpendicular to the Galactic plane. 
The very large error bars for IGR J16479-4514 (Fig. 6a) are 
due to the proximity of the bright low-mass X-ray binary 
GX 340-1-0. We therefore cannot draw conclusions about 
the fluxes of the latter two objects outside of outburst. The 
counts shown are normalised to be the expected 2-10 keV 
count rate if all five detectors making up the PCA were 
turned on. In practice, from one to five detectors are on at 
any given time, and this will be reflected in varying error 
bars from point to point. In these units, the Crab would be 
on the order of 10000 counts s~^. 

The second panels (b) of Fig. 1-8 show the Lomb- 
normalised periodogram of the data, used to search for or- 
bital periods. The gray horizontal line represents the largest 
expected value for each period in 10* randomised trials (the 
real data values were randomly reassigned to the real data 
acquisition times in each trial and the periodogram taken). 
The output of the simulations has a slightly different max- 
imum value for each value of period, but there is no signif- 
icant dependence of the sensitivity on period; the variation 
looks like noise and does not maintain repeatable features 
when the initial random number seed of the full set of ran- 
domised trials is changed. Thus we replace the noisy sensi- 
tivity level of the simulations with a constant value equal 
to the average of all the period bins to give the line shown, 
which is the 99.99 per cent significance threshold for the case 
of a known period. 

The third panel (c) of each figure shows an integral 
luminosity histogram (the fraction of measurements above 
each fiux value), with the fluxes calculated from the count 
rates using the PCA instrument response and the distances 
shown in Table 1 (where more than one distance is cited, 
the first is used). Systems with no good distance estimate 
(XTE J1901+014 and XTE J1743-363) have been placed at 
3 kpc, a distance typical of the estimates for sources that 
have them. If the nature of these two objects is completely 
different from that of the supergiant binaries, their distances 
could be quite different as well. For example, if either is 
a cataclysmic variable, the distance and luminosity could 
be much less. Observation of most of these sources in the 
Galactic bulge is probably due to the region being well ob- 
served rather than any physical association with the bulge 
populations; note that the distances derived for the sources 
with known companions are all well in the foreground of the 
Galactic Centre (Table 1). In these panels, the gray lines 
shown for comparison are the integral of a normal distribu- 
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Figure 1. X-ray time series analysis for XTE J1739-302. Top (panel a): X-ray count rate as a function of time from the bulge scan 
programme. The rate shown is the equivalent rate for all 5 detector units (PCUs) summed over the 2-10 keV band, even though fewer 
detectors were often operating, particularly later in the mission. Bottom left (panel b): Lomb normalised periodogram of the data in 
panel a. The grey line represents the average maximum value in each period bin in 10* randomised trials. Bottom middle (panel c): 
histogram of the fraction of samples above a given luminosity as a function of luminosity, using the first-listed distances from Table 1 
(see text). Bottom right (panel d): autocorrelation function of the data in panel a. 




Figure 2. Data products as described in Fig. 1, for IGR J17544-2619. 



tion centred at zero whose standard deviation is equal to 
the average error bar for the lowest 20 per cent of flux ob- 
servations from each source. This represents the expected 
distribution if there were no real fiux from the source but 
only random noise. 

Finally, the last panel (d) of each figure gives the auto- 
correlation function of the data for each system. To generate 
this curve, we interpolated the time history of each system 
on to a 1 dy uniform grid, setting all the values to zero dur- 
ing the long gaps that occur each year when the Sun passes 
too close to the source. The autocorrelation function was 
calculated on the interpolated grid. Since the source sam- 



pling is usually at approximately 3.5 dy, these plots are only 
informative at times longer than that. 

In Table 2 we present three different measures of the 
percentage of time that each source is active. Considering 
the range of behaviours and luminosities among the sources, 
no single measurement of what it means for a system to be 
"active" or "in outburst" is universally sufficient, nor do we 
expect one of these three definitions to be useful for every 
purpose. The first measure is the percentage of time spent 
with a flux of at least three standard deviations above zero; 
the standard deviation in this case is defined by averaging 
the lowest 20 per cent of the error bars from all the bulge 
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Figure 3. Data products as described in Fig. 1, for 4U 1700-377. 




Figure 4. Data products as described in Fig. 1, for SAX J1818. 6-1703. 



scan fits. These errors are partly statistical and partly based 
on fitting uncertainties from the influence of other nearby 
sources. In general, the lowest 20 per cent of the errors 
correspond to times when the source is undetectable to the 
PCA, with the exception of the sources that are almost 
always visible, such as 4U 1700-377. The second measure is 
the percentage time that the source is above 40 counts s~^ 
(2-10 keV, five detectors equivalent), giving a constant 
threshold in observed flux. The last measure is the percent- 
age of time spent at a luminosity above 1x10"^^ erg s~^ at 
the first listed distance to each source shown in Table 1. 
Bulge scan count rates (2-10 keV) are converted to energy 
fiux (also 2-10 keV) by assuming a thermal bremsstrahlung 
spectrum with temperature 20 keV, typical of SFXTs, 
and using the WebPIMMS tool at NASA's High Energy 
Astrophysics Science Archive Research Centre (HEASARC) 
I http://heasarc.nasa.gov/Tools/w3pimms.html retrieved 



18 April 2011). Any of these can be used as a simplified 
activity parameter, but more information is contained in 
the full luminosity histograms (panel c of Fig. 1-8). 

Table 2 also gives the average and maximum luminos- 
ity of each source, the average value for the bottom half of 
luminosity measurements, and the average of the bottom 
20 per cent of luminosity uncertainties. The latter number 
serves two purposes: first, to make it clear when the average 
of the bottom half of luminosity is significantly above zero 
(only for 4U 1700-377 and XTE J1901+014), and second, to 
give the luminosity threshold used for the first activity in- 
dex: e.g., XTE J1739-302 spent 4.8 per cent of its time with 
a fiux above 3 x (2.30 x 10^^ erg s~^). The negative aver- 
age luminosities for the bottom half of observations in IGR 
J17544-2619 and IGR J16479-4514 are due to the statistical 
and systematic uncertainties in the background-subtracted 
PCA count rates. Negative values in these cases just indi- 
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Figure 5. Data products as described in Fig. 1, for AX J1845. 0-0433. 
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Figure 6. Data products as described in Fig. 1, for IGR J16479-4514. 



cate that in a typical measurement, the uncertainties due to 
subtracting both instrumental background and the influence 
of other sources are greater than the typical count rate of 
the source. 

iRomano et al.l l|201lh showed luminosity histograms for 
three sources (diflerential instead of integral), but sum- 
marised these results with an inactivity measurement in- 
stead. That was most appropriate for their data set, which 
is only 2 years long - and therefore contains fewer outbursts 
- but very sensitive, since the Swift XRT is a true imaging 
instrument. Thus they have good discrimination between 
low and very low emission levels. The RXTE PCA bulge 
scan data are not as easy to interpret at the lowest fluxes, 
even though PCA is also very sensitive, due to complica- 
tions in interpreting the scan data due to the structure of 
diffuse X-ray emission in the Galactic ridge. On the other 
hand, the PCA data set contains over 4000 pointings cov- 



ering over 12 yr for some of the sources we include, so that 
our database is the most sensitive to the distribution of the 
rarer high fluxes. 

3 RESULTS AND DISCUSSION 
3.1 XTE J1739-302 

iDrave et al.l (|2010l ) recently reported evidence for a 51.47 d 
period in this system, based on two separate analyses of IN- 
TEGRAL/IBIS data. In the first, they performed a Lomb- 
Scargle periodogram on the full IBIS light curve and found 
the highest power at this period, with a significance greatly 
exceeding 99.999 per cent confidence. They also collected 
a number of outbursts from IBIS, Swift, and ASCA data 
and plotted their frequency of occurrence and intensity with 
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Figure 7. Data products as described in Fig. 1, for XTE J1743-363. 



XTE J1901+014 
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Figure 8. Data products as described in Fig. 1, for XTE J1901+014. 



Table 1. System characteristics for eight fast transients 



Name 


Companion 


Distance 


Period 


XTE J 1739-302 


08Iab(f) (1,2) 


2.7 (2) 


51.47 ± 0.02 (3) 


IGR J17544-2619 


09Ib (4) 


3.6 (2,4) 


4.926 ± 0.001 (5) 


4U 1700-377 


06.5Iaf+ (6) 


1.9, 2.1 (6,7) 


3.41161 ± 0.00005 (8) 


SAX J1818.6-1703 


~BOI (9,10) 


2, 2.1 (10,11) 


30.0 ± 0.2 (12,13) 


AX J1845.0-0433 


09Ia (14) 


7 (14) 


unknown 


IGR J16479-4514 


08.5I, 09.5Iab (2,15) 


4.9, 2.8 (2,15) 


3.3194 ± 0.001 (16,17) 


XTE J 1743-363 


unknown 


unknown 


unknown 


XTE J1901+014 


faint in K (10) 


1-7 kpc (10,18) 


unknown 



(1) Negueruela et al. 2006; (2) Rahoui et al. 2008; (3) Drave et al. 2010; (4) Pellizza, Chaty, & Negueruela 2006; (5) Clark et al. 2009; 
(6) Ankay et al. 2001; (7) Megier et al. 2009; (8) Hong & Hailey 2004; (9) Negueruela & Schurch 2007; (10) Torrejon et al. 2010; (11) 
Negueruela, Torrejon, & Reig 2008; (12) Zurita Heras 8z Chaty 2009; (13) Bird et al. 2009; (14) Negueruela et al. 2007; (15) Nespoli, 
Fabregat, &; Mennickent 2008; (16) Jain, Paul, & Dutta 2009; (17) Romano al. 2009; (18) Karasev, Lutovinov & Burenin 2008 
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Table 2. X-ray luminosity parameters for eight fast transients 



Name 


Per cent 
> 3cr 


Per cent 
> 40 c/s 


Per cent 

> 1 X 10'-^'' erg s-l 


Average" 


Average, lower 
50 per cent" 


rms, lower 
20 per cent" 


Maximum" 


XTE J 1739-302 


4.8 


4.1 


2.3 


2.25 


-0.15 


2.30 


115.86 


IGR J17544-2619 


5.8 


2.6 


3.6 


1.25 


-2.10 


2.41 


156.58 


4U 1700-377 


85.3 


84.4 


77.1 


86.79 


15.91 


1.16 


1240.55 


SAX J1818. 6-1703 


26.4 


15.2 


3.1 


2.58 


0.75 


0.96 


39.40 


AX J1845.0-0433 


35.5 


6.5 


68.4 


27.00 


5.56 


7.10 


384.19 


IGR J16479-4514 


11.2 


42.3 


56.6 


22.80 


-9.21 


23.59 


193.66 


XTE J 1743-363 


50.6 


10.0 


31.4 


10.36 


1.97 


1.97 


297.95 


XTE J1901+014 


83.9 


1.9 


0.6 


5.20 


3.85 


1.13 


10.90 



" luminosity in units of lO"^* erg s ^ at the first listed source distance in Table 1, or 3 kpc for XTE J1743-363 and XTE J1901+014- 



respect to orbital phase, finding that most outbursts hap- 
pened around phase 0.5 according to their ephemeris with 
zero phase defined arbitrarily at MJD 52698.2. 

We find no evidence for an y power at this pe riod in 
our periodogram (Fig. lb), but iDrave et all (|2010i' ) point 
out that this was also true of Swift and RXTE ASM data. 
We note that the PCA bulge scan data presented here have 
a longer baseline than either Swtft or INTEGRAL and a 
much higher sensitivity than the RXTE ASM. The twice- 
weekly monitoring means that the 51.47 d period is very 
well sampled throughout the 11 years. 

Our result does not, howeve r , nece ssarily constitute a 
contradiction of the iDrave et al.1 (|2010h result. Each PCA 
measurement is a snapshot lasting only a few seconds as the 
spacecraft scans ac ross the source. Mo st INTEGRAL mea- 
surements used by iDrave et al.1 l|2010t ). on the other hand, 
are full 2000 s science windows of that instrument. If we 
hypothesise that the typical behaviour of the source near 
periastron is to produce small fiares that tend to come at 
least once every 2000 s but with a duty cycle of << 50 per 
cent, then it would be plausible that the PCA should miss 
making positive detections at many periastron transits while 
INTEGRAL might catch most of them. 

We the r efore turn to the other technique used by 
iDrave et al.l (|2010l ). the identification of the phase of 
medium-to-large outbursts according to the ephemeris they 
deduced from their periodogram. We have identified 48 
outbursts th at were not included in the sample given by 
iDrave et al.l (2010). These include events from the bulge 
scans presen t ed he re, including some that were presented in 
ISmith et al.l (|2006l ) . outbursts observed during PCA point- 
ings to IE 1 740.7-2942 given in th e same work, Smft out- 
bursts from iRomano et al.1 (l2009l ). and one from INTE- 
GRAL /IBIS that was not included bv iDrave et al.l l|2010l ). 
We also include the m ajor outburst observed by ASCA in 
2002 (MJD 51248.3) jSakano et al] l2q02h. which was in - 
correctly assigned a phase of 0.468 bv iDrave etHI (|2010l ). 
near the apparent periastron, but was not used in the his- 
togram in their fig. 4, which included 35 events from INTE- 
GRAL /IBIS alone. The correct phase of the ASCA outburst 
is 0.83. In many cases, more than one spacecraft observed 
what was probably the same outburst. Since there is no 
prescribed way to distinguish one long outburst from two 
separate ones, we somewhat arbitrarily discarded all new 
outbursts that were within 1 dy of one of those used by 
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Figure 9. Histograms of out burst phases for X TE J1739-302 
according t o the ephemeris of lOrave et al. I 1 I2OIOI) . A; outbursts 
reported bv lDrave et al. I bold ). B: new sample (see text and Ta- 
ble 3). C: the sum of panels A and B. D: the brightest ten events 
in the new sample. 



iDrave et al.l (|201G| ). and similarly made sure that no two of 
the new data points were within 1 dy of each other. The new 
events are given in Table 3. 

Fig. 9 shows histograms o f outburst phase fo r the orig- 
inal sample of events used bv iDrave et al.l (|2010| ). the new 
events, and the sum. The h istogram in Fig. 9 a is identi- 
cal to that given in fig. 4 of [d rave et all l|2010l ). The new 
sample, in Fig. 9b, also peaks near 0.5 in phase. To compare 
these results to the null hypothesis, we calculated the sum of 
the squares of the deviations from the average for each plot 
and compared it to the same quantity for 10® simulations 
of the null hypothesis with the same number of events. The 
real data sets from panels A, B, and C had a sum greater 
than or equal to all but 5.6, 40.4, and 2.4 per cent, respec- 
tively, of the corresponding simulations. This shows that 
while the new data set provides no evidence for periodic- 
ity in its own right, it does increase, and not decrease, the 
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Table 3. Outbu rsts 
iDrave et aEI | |2010D 



of XTE J1739-302 not included by 



MJD 


Pll3,S6 


A/I&jor Outburst? 


Source 


50672.7 


0.647 


Y 


RXTE pointed 


51248.3 


0.830 


Y 


ASCA 


51849.2 


0.505 


Y 


RXTE pointed 


51977.9 


0.005 


Y 


RXTE scan 


52452.9 


0.234 


Y 


RXTE pointed 


52888.3 


0.693 


Y 


INTEGRAL/lBlS 


53569.7 


0.932 


Y 


RXTE scan 


51265.8 


0.170 




" 


51338.5 


0.582 






51660.8 


0.844 






51677.7 


0.173 






52107.5 


0.524 






52560.3 


0.321 






52779.7 


0.583 






52945.1 


0.796 






53140.5 


0.594 






53301.8 


0.727 






53646.1 
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significance of the original result when the samples are com- 
bined. Fig. 9d shows the phase histogram of the ten bright- 
est outbursts in the new sample, corresponding to fluxes 
over 6x10"^*^ erg s~^cm~^ or 5xl0^^erg s~^ at 2.7 kpc and 
marked "Major" in Table 3. There is no apparent clustering 
in phase for these events, which include the ASCA outburst 
at phase 0.83, near the minimum occurrence phase for most 
outbursts. 

The integral luminosity distribution at high fluxes 
(Fig. Ic) has a more or less power-law-like distribution (in- 
dex of approximately -1.5), a characteristic shared with 



IGR J17544-2619 (Fig. 2c), the other "canonical" SFXT 
system, and SAX J1818. 6-1703. The autocorrelation func- 
tion (Fig. Id) shows no significant signal. If outbursts occur 
randomly and last for less than the average sampling period 
(3.5 dy), this is to be expected. 

3.2 IGR J17544-2619 

The orbital period reported bv lClark et"aLl l|2009l ) at 4.926± 
0.001 d appears clearly in the PCA data (Fig. 2b) at the 
same period and with comparable error. This signature 
also appears to a lesser extent in the autocorrelation data 
(Fig. 2d), where it might be taken for noise if the period were 
not already known. The intensity distribution and lack of 
significant autocorrelation of nearby observations look very 
similar to XTE J1739-302. If there is indeed an order of mag- 
nitude difference between the orbital periods of XTE J1739- 
302 and IGR J17544-2619, it is rather extraordinary that 
their behaviours are so similar, in frequ ency of outburst s 
(shown here), frequency of inactivity (Romano et al.|[201ll '). 
and in the detailed app earance of the outbursts them- 
selves jSidoli et al.l |2009|. e.g.). The optical co unterparts 
of IGR J17544-2619 JPellizza. Chatv fc Negueru ela 200^ 
and XTE J 1739-302 (|Negueruela et al."2006l: iRahoui et all 
are very similar as well, which argues against an ex- 
pectation that large differences in the structure and in- 
tensity of their winds can compensate for the different 
orbital radii, producing similar X-ray behaviour. Models 
that rely on unusual magnetic properties of the neutron 
star to explain SFXT outbursts (iGrebenev fc Sunvaev|[2007l : 
IBozzo. Falanga fc Stellall2008l : iLi fc Zhangj|201ll ) should be 
studied further to see if they can produce uniform variability 
behaviour under greatly differing wind conditions. 



3.3 4U 1700-377 

This persistently bright wind accretor is shown as a point 
of comparison for the more transient systems. Its maximum 
luminosity is higher than that of the two canonical SFXTs 
discussed above by factors of 11 and 8, but the average lu- 
minosity is higher by factors of 39 and 69. The comparison 
of the average luminosity is probably very good, but com- 
paring the maximum luminosity is harder, since we don't 
know for how long each system tends to stay at this maxi- 
mum in each outburst, or how common they are. The true 
maximum level within this 11-year period may have been 
missed in one system and caught in another. For example, 
we know fr om coverage of one o utburst of IGR J17544-2619 
by Suzaku (|Rampv et al. I [20091 ) that the system reached at 
least 5.3 x 10^® erg s~^. This is as bright as all but one of 
the data points from the bulge scans of 4U 1700-377. The 
fast transients are capable of reaching almost the maximum 
luminosity of this bright system, but spend most of their 
time at low luminosities. 

The well-known orbital period of 4U 1700-377 (see Ta- 
ble 1) is clearly visible, and the autocorrelation function 
shows coherence on a time-scale of ~ 30 dy. This could be 
due to a process intrinsic to the companion star producing 
slow variations in its wind, a feedback between the accre- 
tion luminosity of the neutron star and mass loss from the 
companion, or the presence of an accretion disc around the 
neutron star, with a viscous time-scale of this order. 
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3.4 SAX J1818.6-1703 

This system (|in't Zand et al.lll998h looks similar to the two 
canonical SFXTs in terms of the power-law shape of the 
distribution of luminosities (Fig. 4c). The orbital period 
(See Table 1) and its harmonics are clearly seen in the pe- 
riodogram; even the autocorrelation function shows a cor- 
relation at the orbital period, which we believe represents 
the high probabilit y of an unre l ated outburst at the next 
periastron passage. ISidoli et al.l ()2009l ') found that the out- 
bursts of this system, studied in depth with Swift, show very 
similar characteristics to the canonical SFXTs. Its 30 dy 
orbital period is approximately halfway between those of 
IGR J17544-2619 and XTE J1739-302, showing again that 
very similar outburst duty cycles and luminosities are seen 
with different orbital parameters. While the system has been 
shown to have as d eep a quiescence as the canonical SFXTs 
IBozzo et al.l l20o3). there has not been enough monitoring 
by a sensitive telescope such as Swift to get the good data 
on frequ ency of quiescence th at are now available for other 
systems omano et al.|[201ll ). 

3.5 AX J1845. 0-0433 and IGR J16479-4514 

There was only one major and one minor outburst of 
AX J1845. 0-0433 caught by the bulge scans, as can be seen 
in Fig. 5, panel (a), and no obvious outburst of IGR J16479- 
4514. The contamination of AX J1845.0-0433 by Galactic 
plane diffuse emission, visible as a bifurcation between those 
scans taken perpendicular and parallel to the plane, makes 
any conclusion about the flux outside of outbursts difficult; 
the contamination of IGR J16479-4514 by GX 340-1-0 pro- 
duces a similar problem. While both of these sources show a 
fairly high average luminosity in the fifth column of Table 2, 
we do not claim this result as significant; repeated obser- 
vations with a focusing instrument such as the Swift XRT 
would be preferable to establish the quiescent and low-flux 
behaviour of these sources. 



3.6 XTE J1743-363 

3.6.1 Position and identification 

This source was discovered by 

iMarkwardt. Swank fc Marshalll l|l999f ). Its occasional 
fast outbursts have led to the suggestion that it might be 
a SFXT, and the X-ray spectrum in outburst could be fit 
with a thermal bremsstrahlung spectrum with temperature 
22 keV, typical of X-ray pulsars and SFXTs (jSeuera et al.l 
l2006l ). iRatti et al.l (|2010l ). using a Chandra observation of 
the field (2008 Feb 8, principal investigator M. Mendez), 
identified a cluster of 11 X-ray photons corresponding to 
the position of the bright infrared source 2MASS 17430133- 
3622221. This is the only significant source in the Chandra 
field, ahhough it is possible that XTE J1743-363 was 
entirely quiescent at that time and the X-rays come from 
another system. Fig. 10 compiles some X-ray position mea- 
surements for XTE J 1743- 363: the initial RXTE estimate 
at th e source's discovery (jMarkwardt. Swank fc Marshalll 
[liii) (at 10', its error circle is too large to appear on 
the plot), the most precise INTECRAL/IQI S position 
(0.7' error circle, 90 per cent confidence) (|Bird et al.l 



-36.40 





1 1 1 1 1 1 1 1 








OLD\_ 




X I 
/ 1 




\ 

\ 




/ 1 


□ \ 


1 




/ \ 


\ NEW V 


1 




/ \ 


\ RXTBA 


1 












Swift BAT^ 


\ 


/ 






OChcwdra ^ 








X INTEGRAL 






, , , , 1 , , , 





265.85 265.80 265.75 265.70 265.65 

RA 

Figure 10. Summary of position measurements for XTE J1743- 
363. The Chandra position has too small an error circle to be 
plotted on this scale, and the first RXTE position has too large 
a circle. The new RXTE position is shown with the 3(t contour 
(dashed line) as well as the la contour (solid line). See text. 

l2010t ). t he Swift BAT pos i tion from the 58-month cat- 
alogue teaumgartner et al.l |2010| ). the Chandra position 
(jRatti et al.ll2010l ). and our new, refined RXTE position. 
Earlier Swift/3AT positions, while different, are consistent 
with both the newe s t BAT position and th e Chandra source 
dTueller et al.ll2010l : ICusumano et al.ll2010l ). 

The RXTE position was derived from a dedicated scan 
of the source using the PGA instrument. The observation 
was on 1999 Feb 15 at 01:15 UT for 2.4 ksec. The scan pat- 
tern consisted of a connected path, including scans along 
constant right ascension and declination. The model in- 
cluded XTE J 1743-363 and one other known source in the 
field of view, 4U 1746-37, at a fixed position. The primary 
change in deriving the refined position was that the inten- 
sity of XTE J 1743-363 was allowed to vary with a long-term 
linear trend. The uncertainties of individual measurements 
were increased in order to de- weight the strong variability of 
the source, which was near 35 per cent rms. While the Chan- 
dra position is slightly outside the 3a RXTE contour, the 
susceptibility of the RXTE fit to variability, combined with 
the agreement of INTECRAL/IBIS with Chandra (spanning 
both the low- and high- energy ends of the PGA energy 
range), leads us to conclude that all the observatories are 
seeing the same source. 

3.6.2 X-ray characteristics 

In hopes of getting further insight into the nature of this 
source, we examined twelve archived RXTE/PCA pointings, 
some from 2000, when it was often bright, and some from 
2007, when it was much less active (see Fig. 7). Table 4 shows 
the dates, RXTE observation identification numbers, and 
count rates for these pointings. Fig. 11 shows some sample 
light curves and Fig. 12 shows the power spectrum made by 
averaging the power spectra from five of the observations. 
Four observations made on MJD 54115 showed low and 
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Table 4. List of PCA pointings to XTE J1743— 363 and absorbed bromsstrahlung fits 



Observation ID 


MJU 




Kate 


nrl 


K i (ke V ) 


Normalisation'^ 


40408-01-01-00 


51224 


051 


12 


11.05 ± 0.22 


13.76 ± 0.35 


5.509 ± 0.082 


40408-01-02-00 


51230 


781 


22 


6.11 ± 0.16 


13.67 ± 0.34 


8.24 ± 0.11 


50138-02-01-00 


51706 


236 


23 


17.40 ± 0.43 


14.66 ± 0.54 


11.44 ± 0.27 


50138-02-02-00 


51707 


631 


13 


10.55 ± 0.51 


15.74 ± 0.97 


4.64 ± 0.15 


50138-02-03-00 


51710 


759 


11 


10.26 ± 0.66 


11.42 ± 0.77 


4.33 ± 0.20 


50138-02-04-00 


51713 


016 


11 


8.77 ± 0.53 


10.54 ± 0.58 


4.01 ± 0.16 


92047-04-01-02 


54115 


404 


2.0 








92047-04-01-04 


54115 


461 


2.0 








92047-04-01-00 


54115 


520 


1.9 








92047-04-01-01 


54115 


910 


1.9 








92047-04-02-00 


54301 


307 


5.3 


11.6 ± 1.4 


14.5 ± 2.2 


1.57 ± 0.14 


92047-04-02-01 


54303 


272 


3.2 


7.4 ± 1.7 


15.9 ± 3.8 


0.503 ± 0.076 


92047-04-02-02 


54304 


383 


3.2 


8.4 ± 1.6 


16.4 ± 3.7 


0.550 ± 0.075 



" counts/s/PCU. 
Xl022 cm-2. 
xlO~^ as defined in equation (1). 



Table 5. Absorbed bremsstrahlung plus blackbody fits to XTE J 1743-363 spectra 



Observation ID nH" Brcmss. kT (keV) Brcmss. norm.* Blackbody kT (keV) Blackbody norm.'' 

40408-01-01-00 8.52 ± 0.47 29.5 ± 6.7 3.09 ± 0.30 1.673 ± 0.042 0.79 ± 0.23 

40408-01-02-00 3.69 ± 0.42 23.5 ± 4.5 4.67 ± 0.50 1.801 ± 0.050 1.21 ± 0.33 

50138-02-01-00 13.44 ± 0.86 24.9 ± 6.2 6.15 ± 0.89 2.058 ± 0.085 1.82 ± 0.54 



Xl022 cm-2. 

xlO"'^ as defined in equation (1). 

xlO"'' normalised to 10^^ erg s"! at 10 kpc jArnaud et al.ll201lf) . 



50 
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Figure 11. X-ray variability at three representative luminosi- 
ties of XTE J1743-363. Top trace: observation 50138-02-01-00. 
Middle trace: observation 92047-04-02-00. Bottom trace (grey): 
observation 92047-04-01-00. See Table 4. 




0.0001 1 . . . . 

0.001 0.010 0.100 

Frequency (Hz) 



Figure 12. Power spectrum combining the data from observa- 
tions 40408-01-02-00, 50138-02-01-00, 50138-02-02-00, 50138-02- 
03-00, and 50138-02-04-00. No pulse period is visible between 2 s 
and 1024 s. 



nearly identical count rates of 1.9-2.0 counts s~^ and no 
variability. The flux varied during all the other observations 
on time-scales of seconds to minutes. We used xspec version 
12.5.1 to fit the PCA spectra. These four quiet observations 
included an iron line of high equivalent width, and no de- 
tectable interstellar absorption above 3 keV, in contrast to 
all the other observations, which showed significant and vari- 
able absorption. Finally, the count rates of the four spectra 



from MJD 54115 are comparable to the flux level that is 
identified as diffuse from the PCA bulge scans in this region. 
Taking all this evidence together, we believe it quite likely 
that during these four observations we were seeing only the 
Galactic diffuse component and not the source itself. The 
best-fitting continuum spectrum to these four observations 
together was an optically thin thermal bremsstrahlung spec- 
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trum of temperature 10.10 ± 0.50 keV and normalisation 
(5.99 ± 0.16) X 10-^ defined as 



3.02 X 10" 

47rD2 



neiiidV 



(1) 



where D is the distance to the source in centimetres and rie 
and ni axe the el ectron and ion dens ities per cubic centime- 
tre at the source (jArnaud et al.l[2"oilf l. An iron line of equiv- 
alent width 0.43 keV and centre energy (6.63±0.54) keV was 
also required (at 3.7a), and was consistent with being nar- 
row, although the width was poorly constrained (<2.9 keV). 
In fitting the rest of the spectra, we included these two spec- 
tral components with no free parameters. No further iron 
line was needed in any of the other fits. 

In Table 4 we show the spectral parameters of the 
remaining flux, presumed to be from XTE J1743-363 it- 
self. For most of the spectra, an absorbed optically thin 
thermal bremsstrahlung spectrum provided a good fit. We 
used the Tuebingen- Boulder ISM absorption model (tbabs) 
I Wilms. Allen, fc McC"ravll2000l) . For the three spectra with 
the best counting statistics, however, an additional compo- 
nent was necessary to obtain x-square values consistent with 
chance. A thermal blackbody spectrum served this purpose 
well for all three spectra, as it typically does for accreting 
neutron stars. Table 5 shows the resulting fit parameters in 
this case. Including the blackbody component has the effect 
of raising the temperature of the bremsstrahlung component 
and reducing the best-fitting absorption column. 

It is apparent from Tables 4 and 5 that the absorption 
column is both high and variable, and therefore local to 
the system. This suggests a companion with a dense wind, 
which can be characteristic of both blue supergiants and red 
giants such as AGB stars. The absorption column doesn't 
seem to correlate with luminosity at all, suggesting that the 
absorption is taking place in foreground structures in the 
wind rather than in the concentration of material accreting 
on to the compact object at the moment. 
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Figure 13. 2MASS J-band image of the region around 
2MASS 17430133-3622221. The arrow marks the source. 



3.6.3 Optical spectroscopy 

iRatti et al.l (2010) used the infrared colours of 
2MASS 17430133-3622221 to constrain the spectral 
type of the companion, but not severely; they found 
that a giant of type GO-6, a G or K supergiant, or a 
supergiant earlier than 09 were all consistent with the 
2MASS colours, assuming varying degrees of reddening. 
The 2MASS magnitudes are J=9.616, i/=8.305, and 
K=7.624:. The 2MASS source (Fig. 13) is coincident with 
USNO-Bl.O 0536-0518280, with R=16.4 and 16.77 (Rl and 
R2 measurements) and 7=13.78 (not seen in the USNO B 
band) and with the Midcourse Space Experiment infrared 
source MSX6C G353.3697-03.4195, seen only in band A at 
0.1386 Jy. Band A is centred at 8.28 /xm with half-maximu m 
points at 6.8 and 10.8 fim according to lEgan et al.l (|2003l ). 

We obtained a far-red spectrum of the counterpart us- 
ing the Andalucia Faint Object Spectrograph and Camera 
(ALFOSC) on the 2.6 m Nordic Optical Telescope (NOT) in 
La Palma, Spain, on the night of 2009 June 27. The instru- 
ment was equipped with a thinned 2048 x 2048 pixel E2V 
CCD, covering a field of view of 6.'4 x 6!4 with a spatial scale 
0'.'l9/pixel. We used grism #4 with a I'.'O slit, providing a 



I 



Figure 14. NOT finder image covering the same area as Fig. 13. 
Note that 2MASS 17430133-3622221 is resolved into two sources. 
The brighter source to the left (marked with a circle) is the M 
giant discussed in the text. 



resolving power R « 450. In order to remove fringing, an 
internal fiat was taken at the target position. 

Because of the Southern declination of the target, it 
was observed at a very high airmass {Z = 2.4). However, as 
the finder image (Fig. 14) shows, image quality was good. 
2MASS J17430133-3622221 can be clearly resolved into two 
objects, separated by ~ 2", both of which fell within the slit. 
The fainter object is a normal G/K star. The brighter star 
is a very red object, and the obvious counterpart to the 
2MASS source. Its spectrum (Fig. 15) is dominated by deep 
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Figure 15. ALFOSC spectrum of 2MASS 17430133-3622221/USNO-Bl.O 0536-0518280. The positions of some relevant molecular 
bandheads are indicated. Tlie strongest telluric features are also marltod. 



molecular bands. In particular, we see deep bandheads of 
TiO, typical of M-type stars. 

iNegueruela et al] |20l3) have shown that the depth of 
the TiO 8859A bandhead is a very good indicator of effective 
temperature. Though the resolution of our spectrum is too 
low to use it as a T^g calibrator, the fact that its depth is > 
0.5 of the continuum immediately shortward indicates that it 
is not a dwarf star. Only giants and supergiants present such 
deep molecular features. It also indicates that the spectral 
type is later than M6. This is confirmed by the detection of 
VO band heads, which are only visible at spectral types M7 
and later (|Turnshek et al.ll 19851 : [Crav fc Corballvll2009l ). 

The deep depression of the continuum in the 7400 - 
7550Aand7860-8050A regions can only be due to the pres- 
ence of strong VO bands of the extreme red system, as there 
are no telluric features in these two regions. Moreover, the 
VO bandhead at 8624A is also likely detected. Though a 
better spectrum is needed to obtain an accurate spectral 
type, the data are not very different from the expectation 
for spectral type M8. A s supergiants with su ch a late spec- 
tral type are very rare (jLevesque et al.ll2005l l. the object is 
very likely a giant. 

The fundamental properties of such late stars are not 
well calibrated. For M8III stars. Ivan Belle et all IHii) es- 
timate Taft = 3050 K, though this is an extrapolation from 
earlier spectral types. Their fit to the intrinsic (V^ — _ft')o 
colour indicates values between (V — K^a = —8 . and 
iy - K)o = -8.5 at this temperature. As lThe et akl (|l990l ') 
find Mv = -1-1.6 on average for M8III giants, their infrared 
absolute magnitude is Mk ^ —6.4. Late M-type giants may 
be close to tip of the red giant branch (RGB) or on the 
asymptotic giant branch (AGB). Many giant s later than 
M6 ar e brighter than the tip of the RGB (e.g.. iTabur et ahl 
l2009al ). which has observationally been set at Mk = —6.85 



dTabur et al.ll2009bl ). and so most M8 III giants are likely on 
the AGB. 

Theoretical magnitudes calculated from PHOENIX mod- 
els and calibrated against observational d ata suggest that 
all la te-M giants have (J — K)o ~ 1.3 (|Kucinskas et al.l 
!2005l), in good agreement with extrapolation of the colour 
calibra tion of [Straizvs fc Lazauskaitd (|2009l ). The calibra- 
tion of lGroenewegenl ~ 20041 '). based on older models, gives 



(J — K)o = 1.37 for dust-free M8 stars. This reference also 
gives (J — K)o ~ 1.97 for dust-enshrouded M7 stars. Based 
on this, our target cannot be dust-enshrouded, unless the 
insterstellar extinction is negligible. This possibility is ruled 
out by the non-detection in all MSX bands except for A, and 
so we assume that the star is not surrounded by dust. There- 
fore its colour excess is E{J — K)o ~ 0.6 — 0.7. For a standard 
reddening law, Ak ~ 0.4 — 0.5 and so the distance modu- 
lus should be DAI « 13.5, i.e., d « 5 kpc. As observations 
of LMC giants suggest that Mk may deviate by ±0.6 mag 
from the average at a given (J — K) l|Cioni et al.ll200ll ). we 
may settle for d « 5 ± 1 kpc, though there are many other 
sources of uncertainty not considered in this error estimate. 

For instance, these calculations assume that the star has 
characteristics typical of M-type giants in the solar neigh- 
borhood. Our upper distance limit, however, brings us close 
to the distance to the Galactic bulge. At 6 = — 3?4, the 
extinction to our s ource is not incom patible with Galactic 
bulge membership (|Dutra et al.ll2003l ). even if we take into 
account the possibility of bluer colours for bulge M-giants 
(Frogcl 1988,). However, for a distance of 8.5 kpc, the source 
would have Mxg ~ —7.8, o nly compatible with the very 
brightest bulge M-type stars ([Glass fc Schultheij2002l ). Our 
target would then have to be a thermally-pulsating AGB. 

Even in the most likely case that the star is foreground 
to the bulge, it is likely to be variable, as most late-type 
giants are long-period variables. Their amplitude in the K 
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Figure 16. Data products as described in Fig. 1, for thie symbiotic binary GX 1+4. 



band is generally small and unlikely to dominate our error 
budget. On the other hand, photometric determination of a 
pulsation period would be the most direct way to evaluate 
the stellar luminosity and hence distance. 



troscopy of the counterpart rules out this possibility, unless 
the correct companion is the faint, nondescript G/K star 
nearby. It is unlikely, however, that a star as unusual as the 
M8 giant would coincide with the X-ray source by accident, 
while it is quite likely that the dwarf star would do so. 



3.6.4 XTE ,11743-363 as a symbiotic X-ray binary 

The spectral identification of the companion to XTE J1743- 
363 as a late M giant puts the system in the small group of 
symbiotic X-ray binaries, although we have not seen an addi- 
tional hot component to the optical spectrum, which was the 
origin of the term "symbiotic" . For comparison, we show in 
Fig. 16 the PCA bulge scan data for GX 1-1-4, a known sym- 
biot ic X-ray bina ry in which an X-ray pulsar orbits an M gi- 
ant. [Hinkle[er|al] ( 2006 ') used infrared spectroscopy to derive 
an orbital period of 1161 dy, a luminosity of 1270 Lq, and a 
distance of 4.3 kpc, which we assume to scale the luminosity 
distribution in Fig. 16. Comparing GX 1-1-4 with the data for 
XTE J1743-363 in Fig. 7, we see two things in common: a 
noticeable autocorrelation out to 100 dy (while even the very 
bright blue supergiant system 4U 1700-377 only correlates 
out to 30 dy) and a long quiet period (for GX 1+4, the years 
2007 and 2008 and for XTE J1743-363 the last 5 yr). These 
common characteristics lend support to the spectroscopic 
evidence that XTE J1743-363 is also a symbiotic binary. 
Strong emission lines in the optical spectrum and pulsations 
of GX 1+4 have be en taken as evidence of an accretion disc 
IChakrabartv fc Ro che 1997; Jablonski et al. 1997). While 
these characteristics have not been seen in XTE J1743-363, 
we note that our far-red spectrum of XTE J1743-363 was 
taken in the middle of the extended X-ray quiet period still 
in progress. It is possible that a disc forms when the source 
is X-ray active. 

The time variability of XTE J1743-363 also resem- 
bles that of a member of another class of binary entirely: 
the ultracompact binary 4U 0513-40 in globular cluster 
NGC 1851. It also shows highly variable emission, corre- 
lations over many days, and a relatively quiet p eriod of ap- 



3.7 XTE J1901+014 

This system showed a short (< 3.15 hr), bright 
(0.9 Crab) X-ray outburst a t the time of its dis- 
covery by iRemillard fc SmithI (|2002l'). It has since 
shown further outbursts I Krimm et al.l |2010| ) but 
seems to have a relatively high flux in quies- 
cence as well (iKarasev. Lutovinov fc GrebenevI l2007l : 



q ui 

proximately three years (|Maccarone et aUl201Cl ). Our spec- 



iKarasev. Lutovinov fc BureninI l2008f) compare d to the 
SFXTs (|Smith et all l2007l : iKrimm et al.l boid l. It does 
not have an IR/optical counterpa rt bright enough 
to be a supergiant in our Galax y (|Smith et al.l l2007l : 
IKarasev. Lutovinov fc BureninI [2OO8I ) . Although we didn't 
see an outburst in the bulge scan data, we confirm the 
presence of a nonzero quiescent baseline that is not con- 
sistent with the behaviour of t h e SFX Ts in our sample. 
IKarasev. Lutovinov fc GrebenevI ||2007| ). using pointed ob- 
servations of the RXTE PCA, showed that there is a great 
deal of variability (by more than a factor of two over tens of 
seconds) even outside of outbursts. They quoted a typical 
flux of 2.7 mCrab outside of outburst. Our average luminos- 
ity of 5.72 xlO'^* erg s~^ at an arbitrary distance of 3 kpc 
corresponds to 5.3 xlO~^'^ erg s~^ or a bit over 2 mCrab , 
in agreement with IKarasev. Lutovinov fc GrebenevI (|20o3). 
The persistent emission outside of outburst supports the 
more definitive evidence of the lack of a supergiant com- 
panion in separating this source from the class of SFXTs. 
The nature of the system remains a mystery. 
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